Stop-and-go kinetics in amyloid fibrillation 
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Many human diseases are associated with protein aggregation and fibrillation. We present ex- 
periments on in vitro glucagon fibrillation using total internal reflection fluorescence microscopy, 
providing real-time measurements of single-fibril growth. We find that amyloid fibrils grow in an 
intermittent fashion, with periods of growth followed by long pauses. The observed exponential 
distributions of stop and growth times support a Markovian model, in which fibrils shift between 
the two states with specific rates. Remarkably, the probability of being in the growing (stopping) 
state is very close to 1/4 (3/4) in all experiments, even if the rates vary considerably. This finding 
suggests the presence of 4 independent conformations of the fibril tip; we discuss this possibility in 
terms of the existing structural knowledge. 

PACS numbers: 87.14.em 87.15.bk 82.39.-k 
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Protein fibrillation is the process by which misfolded 
proteins tend to form large linear aggregates [l[ . Its im- 
portance is related to the role played in many degenera- 
tive diseases, such as Parkinson's, Alzheimer's, Hunting- 
ton and prion diseases While our knowledge of the 
structural properties of these fibrils improves at great 
pace 0, 0, lH, the dynamics of their growth process is 
still poorly understood. The formation of amyloid fibrils 
involves at least two steps: the formation of growth cen- 
ters by primary nucleation, which is often a slow process, 
followed by elongation through addition of monomers Q . 
In many cases, a so-called secondary nucleation mecha- 
nism is also involved, whereby new grow th centers are 
formed from existing fibrils @, 1, BM EH- Whereas 



the process of secondary nucleation is known to entail a 
number of different mechanisms [TH |. the primary elon- 
gation process has not been elucidated to the same level 
of detail. 

In this letter, we present an experimental and theoret- 
ical study of the elongation process of glucacon fibrils. 
Glucagon is a small peptide hormone consisting of only 
29 amino acids produced in the pancreas. It has the 
opposite effect to that of insulin and therefore increases 
blood glucose levels when released. As a model system 
for protein fibrillation, glucagon kinetics has provided in- 
sights into the early oligomerization stages of the process 
12 . 13, 14, TiL the interplay between growth and fibril 
morphology [Ty, [13 & n( i amyloid branching Here, 
we focus on the properties of the late-stage elongation 
process. 

Experiments were performed on samples of glucagon 



monomers in solution. In order to detect the growth 
a specialized fluorescence microscopy technique was ap- 
plied, the so-called Total Internal Reflection Fluorescence 
Microscopy (TIRFM). This technique utilizes total in- 
ternal reflection to create an evanescent electromagnetic 
field adjacent to the glass slide, thereby exclusively ex- 
citing fluorophores in only a very thin volume. The pen- 
etration depth, d, depends in a specific manner on the 
wavelength and the angle of the incident light, as well 
as the refractive indices of the media [lij . In the setup, 
an argon laser was used along with a fused silica slide 
in contact with water, leading to a penetration depth 
of d = 150 nm. The TIRFM images of the fibrilla- 
tion process were obtained at initial glucagon concen- 
tration of p = 0.25 mg/ml in aqueous buffer (50 mM 
glycine HC1 pH 2.5) with preformed seeds. Images of the 
growth are shown in Fig. [1] at three consecutive times 
t = 0, 216 , 407 min. (see |ll| for details on the experi- 
ment). 

Because the fibrils grow along the glass slide we are 
able to track each fibril length as function of time. We 
monitor 16 independent fibrils for each image frame in 
the experiment. The time interval, At, between frames 
varies from a minimum value of Ai m .; n = f min. to a 
maximum value At max — 35 min., with a typical value 
of At = 10 min.. The total duration of the experiment 
is t = 525 min. The combined results for the 16 fibrils 
are shown in Fig. [5] 

A striking feature of the fibril dynamics is its discrete 
nature, where long periods of growth are interrupted by 
extended periods of stasis (stop state). This prompted 
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FIG. 1: TIRFM images of glucacon fibril growth with initial 
glucagon concentration of p — 0.25 mg/ml in aqueous buffer 
(50 mM glycine HC1 pH 2.5) at three consecutive times after 
the initiation of the aggregation. Red lines in the last picture 
mark examples of fibrils which are tracked during the growth 
process. 



us to collect the statistics of time spent in the growth 
(g) and the stop (s) state, f g (t) and f s (t) respectively, 
for all 16 fibrils. In Fig. these distributions are shown 
on semi- logarithmic plots. Note that the finite sampling 
rate implies that actual time spent in given state can 
only be estimated in terms of upper and lower bounds 
[19( | . The upper estimates are shown by the dashed blue 
curves and the lower estimates are shown by the full blue 
curves. As seen in the figure, the difference between 
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FIG. 2: Length as a function of time for 16 fibrils tracked 
from the images shown in Fig 1. Note the long plateaus, cor- 
responding to the stop states, followed by shorter (on average) 
growing periods. The average growth is indicated with a dot- 
ted blue line. The sampling time between each image can be 
seen from the time separation between each point. 

the distributions for the upper and lower estimates of 
both f g and f s are marginal. All distributions are very 
well fitted by exponential functions, f g (t) ~ exp(— fc_i) 
and f s (t) ~ exp(— fc+i), as shown by the yellow curves 
(here, the dashed yellow curve is the fit to the upper es- 
timates and the full yellow curve is the fit of the lower 
estimates). The fits are of excellent quality over almost 
two decades, as signified by high R values (R 2 > 0.98). 
The values obtained are k + = 9.0 • 10 -3 min -1 and 
k- = 2.8 ■ 10~ 2 min - . A series of four independent 
experiments A — D have been performed with the same 
glucagon monomer concentration (0.25 mg/ml) and pH 
(2.5) with differences in the seed concentrations (~20% 
variation) and data acquisition only. Specifically, the 
maximal frame length was varied by a factor 3 and the 
minimum frame length by a factor 18. In all four cases 
the same analysis and procedure was performed resulting 
in exponential distributions of similar high quality. The 
results are summarized in Table |H 

The observed stop and go behaviour of fibril dynam- 
ics is clearly not associated with the discrete nature of 
monomer attachment (20|. The simplest model of the 
process is to assume that the fibril exhibits two internal 
states: one in which it is allowed to grow, with a rate g 
and one in which it cannot grow. The intrinsic transi- 
tion rates between the two states is then identified with 
the observed transition rates k + (stop^ growth) and fc_ 
(growth — > stop), see Fig. |4j 

The Markovian nature of the model implies, that the 
time spent in each state is exponentially distributed 
which is indeed consistent with the results from the data- 
analysis, Fig. ([3]), leading to estimates of k + and fc_. De- 
noting the total rate k = fc + + fc_ , one can calculate the 
probability of being in the growing state, p + = k+ /k and 
in the stopped state, p_ = k^/k. Table U presents a sum- 
mary of the parameters in 4 different experiments. By 
comparing the four different series, we observe that the 
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FIG. 3: Growth f g (t) (a) and stop f s (t) (b) times distribu- 
tions on semi- logarithmic scales. Blue lines are data and yel- 
low straight lines are the exponential fit. Both for data and 
fits, continuous lines are the lower estimates and dashed lines 
are the upper estimates (see text). All fits are of extremely 
good quality as indicated by the large R-value, R 2 > 0.98. 
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FIG. 4: The simplest two-state model of fibril growth. The 
fibril is assumed to be in one of two states with intrinsic tran- 
sition rates of k+ and In the growing state, monomer 
attachment occurs with the rate g. 



rates fc+ and fc_ vary quite significantly, up to a factor 
3. A possible explanation for this could be the variation 
in sampling frequency (as is indicated in Table [TJ) , due 
to heating of the sample by the laser. Indeed, the de- 
pendency of the fibril growth on both the laser intensity 
as well as the illumination time has been observed un- 
der similar experimental conditions for /^-microglobulin 
kinetics (2lT |. 

The striking result is that, even though both transition 



rates vary among experiments, they combine in such a 
way that the probabilities of growing and stopping, p + 
and p-, do not change appreciably in different exper- 
iments. In particular, p+ is always very close to 1/4 
and consequently p_ is very close to 3/4. Notice that, 
if the difference between the growing and the stopped 
state would have been due to an energy gap, one would 
have expected the population ratio of the two states to be 
much more sensitive to variations in the individual tran- 
sition rates, k±. Conversely, the constancy of this ratio 
suggests that the energy difference between the growing 
and the stopped states to be irrelevant. This ratio could 
then reflect the presence of three stopped configurations 
for each growing one, all of them being isoenergetic. 

One may wonder whether the observed '1/4-3/4' law 
is really independent of the monomer concentration, as 
in the model. Due to experimental limitations, this has 
not been directly tested. However, the observed average 
growth rate, <?, displays more than a two-fold variation 
from Exp. A to Exp. D, suggesting notable differences in 
the local monomer concentration of the growing fibrils. 
If the state probabilities were sensitive to changes in the 
concentration, we would have expected a correlation be- 
tween p± and g. We further tried to test the pertinence 
of the state probabilities by dividing the data sets into 
two parts, one corresponding to the first half of the ex- 
periment and the other to the second half, and measuring 
p + and p_ separately in the two parts. Unfortunately, 
the statistical noise increased significantly when consid- 
ering sub-parts of the data series and neither added to 
nor altered our conclusions. 

To conclude, we have presented a stop-go model for 
glucagon fibril elongation. In fact, similar kinetic behav- 
ior has recently been reported for the fibril elongation 
of Af3 -peptides 



[22|, |23j as well as for a-synuclein 
although the timescales involved are 1-2 orders of mag- 
nitude faster. This suggests that the observed stop-go 
kinetics reflects the presence of some kind of structural 
change at the fibril ends, which is not necessarily spe- 
cific to glucagon. In this picture, the protein properties 
would affect the barrier height, and thus the timescale of 
the process, only. 

If our hypothesis about the existence of approximately 
isoenergetic states holds up to scrutiny, it suggests that 
there may be an additional dimension in the fibrillation 
energy landscape that cannot easily be identified by con- 
ventional techniques. Glucagon is known to adopt a num- 
ber of different conformations depending on the fibrilla- 
tion conditions, but these conformations differ consider- 
ably in energy and are unlikely to co-exist to an equal 
extent pBj . Rather, it is possible that we have a num- 
ber of closely related states with different propagation 
properties which are separated by high local activation 
barriers within a relatively flat ground state level. Fur- 
ther experimental studies are required to establish the 
validity of this suggestion. 
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Parameters 


Exp. A 


Exp. B 


Exp. C 


Exp. D 


k+ [min -1 ] 


0.0090 


0.0192 


0.0075 


0.0059 




±0.0001 


± 0.0004 


± 0.0001 


± 0.0001 


fc_ [min -1 ] 


0.0278 


0.04780 


0.023 


0.0166 




± 0.0005 


± 0.00002 


± 0.002 


± 0.0002 


P+ 


0.244 


0.249 


0.247 


0.260 


P- 


0.756 


0.751 


0.753 


0.740 


g [nm/min] 


135 


116 


111 


56 


Tot. time [min] 


525 


340 


630 


1030 


Max. FL [min] 


35 


27 


24 


90 


Min. FL [min] 


1 


3 


1 


18 


Avr. FL [min] 


13.7 


11.8 


20.7 


30.9 



TABLE I: Estimated kinetic parameters (row 1-5) and pa- 
rameters for the data acquisition (row 6-9) in 4 different ex- 
periments, A-D. FL stands for frame length. All figures of 
this paper refer to the experiment of the first column, Exp. 
A. The average growth rate of the fibrils over the total time 
of each experiment is indicated by g (no appreciable effect of 
depletion is observed). These values are comparable to the 
measurements in flUl. 
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